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We present experimental nitrogen-broadening coefficients derived from Voigt profiles of 
isotropic Raman Q-lines measured in the v 2 band of acetylene (C 2 H 2 ) at 150 K and 298 K, 
and compare them to theoretical values obtained through calculations that were carried 
out specifically for this work. Namely, full classical calculations based on Gordon's 
approach, two kinds of semi-classical calculations based on Robert-Bonamy method as 
well as full quantum dynamical calculations were performed. All the computations 
employed exactly the same ab initio potential energy surface for the C 2 H 2 -N 2 system 
which is, to our knowledge, the most realistic, accurate and up-to-date one. The resulting 
calculated collisional half-widths are in good agreement with the experimental ones only 
for the full classical and quantum dynamical methods. In addition, we have performed 
similar calculations for IR absorption lines and compared the results to bibliographic 
values. Results obtained with the full classical method are again in good agreement with 
the available room temperature experimental data. The quantum dynamical close- 
coupling calculations are too time consuming to provide a complete set of values and 
therefore have been performed only for the R(0) line of C 2 H 2 . The broadening coefficient 
obtained for this line at 173 K and 297 I< also compares quite well with the available 
experimental data. The traditional Robert-Bonamy semi-classical formalism, however, 
strongly overestimates the values of half-width for both Q- and R-lines. The refined semi- 
classical Robert-Bonamy method, first proposed for the calculations of pressure broad- 
ening coefficients of isotropic Raman lines, is also used for IR lines. By using this improved 
model that takes into account effects from line coupling, the calculated semi-classical 
widths are significantly reduced and closer to the measured ones. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Many experimental infrared studies have been devoted 
to the nitrogen broadening of acetylene lines of various 
fundamentals or combination bands [1-9]. Table 1 of the 

* Corresponding author. Tel.: +33223236182: fax: 33223235662. recent studies of Predoi-Cross and collaborators [8,9] gives 

E-mail address: Franck.Thibault@univ-rennesi.fr (F. Thibault). an exhaustive summary. In the present work, we complete 
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Table 1 

Pressure broadening coefficients for the R(0) line at room temperature 
(297 + 2 I<), 173 K and 100 K (as a reference for future use, for instance, 
for Titan's atmosphere). 



Room T 

173 K 

100 K 

Calculated 

120.45 (a) 

175.75 (a) 

259.05 (a) 


119.75 (b) 

177.95 (b) 

264.85 (b) 


121.05 (c) 

176.30 (c) 

258.50 (c) 


121.30 (d) 

177.95 (d) 

264.85 (d) 

Experimental 

116.60 (e) 
119.95 (f) 
116.15 (g) 
110.00 (h) 

161.10 (i) 
168.00 (j) 



(a) CC ThAv, (b) CC MThV, (c) CC/CS ThAv, (d) CC/CS MThV, (e) i / 5 band [ 1 ] 
297 K, (f) v,+v 3 band [5] 295 K, (g) v,+v 5 band 297 I< [2], (h) v,+v 3 band 
[ 8 ], (i) V 4 +V 5 band [ 6 ] 173.2 K, (j) v 5 band 173.4 K [3], 


these studies by performing the very first broadening 
measurement of isotropic Raman Q-branch lines in the v 2 
band at 150 K and 298 K. The motivation of this and 
similar works are well known, with the data obtained 
being of interest to (and finding application in) fields as 
diverse as the study of planetary atmospheres -including 
the Earth - combustion processes or pollution control. 

Due to the interest of the C 2 H 2 -N 2 system many 
theoretical calculations of the nitrogen broadening of 
acetylene lines have also been performed [1,3,4,8-11], 
However, these calculations were almost all conducted 
with the semi-classical method of Robert and Bonamy [12] 
using semi-empirical potentials. For instance, in Ref. [8] 
various isotropic parts of this system were tested in order 
to get a good agreement with the experimental data. A few 
calculations were also performed on the basis of the 
energy corrected sudden approximation, for which no 
potential energy surface was used and for which the 
nitrogen molecule was reduced to a pseudo-atom [13,14], 
The very first realistic dynamical calculations performed 
on an ab initio potential energy surface (PES) [15] were 
done in Ref. [9] using the full classical method of Gordon 
[16,17], In the present study, starting with the same PES 
[15], we are extending this effort. We intercompare the full 
classical method (FC) with semi-classical (SC) methods 
[12,18-20], quantum dynamical calculations, available 
experimental results for infrared lines and our results for 
isotropic Raman lines. In addition, this study completes a 
series of papers [21-23] devoted to a comparison of the 
full classical, the close-coupling (CC) and the Robert- 
Bonamy (RB) methods. In that series it has been shown 
that the traditional RB formalism is unable to provide 
accurate pressure broadening coefficients without adjust- 
ing some intermolecular parameters. More recently, Ma 
et al. [20] have shown that by considering effects from the 
line coupling in the standard [12] RB method, one is able 
to obtain a much better agreement with experimental 
values at least for self-perturbed nitrogen isotropic Raman 
lines. The latter study has been completed by a compar- 
ison of the full relaxation matrix as obtained with this 
refined method with quantum dynamical calculations 
[24,25], Meanwhile, this refined version of semi-classical 
RB model has been adapted to the case of IR lines [26], 


The aim of such comparative studies, including full 
classical ones, is to find well founded alternative methods 
to the quantum dynamical calculations, in particular to the 
CC method, which still remain very time consuming even 
for simple systems like N 2 -N 2 or C 2 H 2 -N 2 . 

In the following sections we briefly recall a few points 
about the methods used. Then the experimental set up and 
reduction of data are given. The results of our calculations 
are compared with our experimental Raman data and, 
finally, with bibliographic infrared data. 

2. Theoretical methods 

Binary collisions and the impact approximation for 
lines without line mixing [27] are the shared features of 
the methods summarized below. We disregard any trans- 
lational effects (velocity or speed dependence) on the 
lineshape. Since the PES [15] is only four dimensional in 
the Jacobi coordinates, our calculations are performed 
within the rigid rotor approximation. We are therefore 
unable to provide any lineshift coefficients as well as to 
discuss any vibrational dependence of the collisional line- 
widths. The latter seems to be quite insensitive to the 
different kinds of transitions involving different vibra- 
tional modes [7] but the situation is not so clear regarding, 
for instance, Fig. 6 of [6], Fig. 13 of [7], Fig. 5 of [8] and Fig.4 
of [9], However, such experimental differences may arise 
because of the different lineshape models used. Regard- 
less, our aim is not to provide accurate pressure broad- 
ening coefficients but to compare the methods we have 
used starting with the same PES. 

2.1. Potential energy surface used 

The calculations have been performed on the four- 
dimensional ab initio PES of Thibault et al. [15]. This PES 
has been expanded over bispherical harmonics (see Refs, 
therein [23]): 

V(R.6^,O 2 ,0) = x A tlJ . 2i (6h,# 2 ,<£) (1) 

where (6fi, e 2 , 0) are the three angular Jacobi coordinates 
which describe the relative orientation of the two mono- 
mers. R is the distance between the two monomer centers 
of mass and defines the intermolecular z axis, 0{ 12) defines 
the bending angle of the monomer axis relative to the 
intermolecular axis (the indices “1” and “2” refer to the 
acetylene and nitrogen molecules respectively) and 0 is 
the rotational angle around the collisional z axis oriented 
from “1” to “2”. The values L lt L 2 and L are even because 
the colliding pair is formed of centrosymmetric molecules 
and |L 1 -L 2 | <L< |Li+L 2 |. In Eq. (1) 85 terms have been 
used, up to (L],L 2 ,L)=(8,8,16). The angular functions 
A L ^ L2L (9- l ,O 2 ,0) are defined in the case of two linear 
molecules by a normalized product of spherical harmonics 
for monomers 1 and 2 as 

AL,L 2 L(,Rl,d 2 ,0) 

/2L+1\ 1/2 

= I,(Ltml 2 -m\LO)Y 1 !'fe,,O)Y[ 2 m (e 2 ,0) (2) 

where V™ and Yf™ are the ordinary spherical harmonics, 

is a Clebsch-Gordan coefficient, and |m| < minfi^f^). 
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The radial coefficients V Li l 2 l(R) were extrapolated, using 
C a /R " functions, from 16 to 30 A since the ab initio points 
end at 16 A. 

Such an expansion of the PES is suitable for the 
MOLSCAT quantum dynamical code [28] and for the classi- 
cal calculations. However, in order to be used in conjunction 
with the RB method, in a fixed reference frame (Fig. 5 of 
[12]), a transformation has to be performed, leading to 

-m o) 

(3) 


2.2. Pressure broadening cross-sections and coefficients 


Considering C 2 H 2 (index “1”) in mixtures with natural 
N 2 (index “2") composed of a 2:1 ratio of ortho-(oN 2 ) and 
para-(pN 2 ) species associated with even and odd rota- 
tional quantum number j 2 , the broadening coefficient 
r(ji’T) (in cm -1 atm -1 ) at a given temperature T can be 
expressed as weighted sums of partial pressure broad- 
ening coefficients [22,23,29,30]: 

2 i 

/Ol > T) = gZj 2 evenPh /O' I J 2 , T) + 3 Zj 2 oddAkl'Ol ’h> 'O- (4) 

In Eq. (4) pj 2 is the normalized rotational population 
and 

rOi J 2 J) =/<v<r (n) (/j ,h,Ekin)) =f vff (n, 0'i ,h, T), (5) 

where tr m (j\j 2 ,T) is the Boltzmann average [29] over 
kinetic energy of the partial pressure broadening cross- 
sections cr (n> (j 1 J 2 ,E k i n ) (in A 2 ) and f = (D/2irc)(273/T), 
where D is the number density at standard temperature 
and pressure (D= 2.69 10 19 molecule cm -3 ) and c is the 
speed of light (in cms -1 ). v is the mean relative speed 
associated to the thermal speed distribution: 
v = y/&kfTJjrp, where p is the reduced mass of the collid- 
ing pair in u units. Finally, the “exponent” (n) denotes the 
order of the tensor of the interaction between matter 
and light (n=0 for isotropic Raman lines and 1 for 1R 
transitions). 

An approximation often employed for the semi-classical 
method (for technical reasons) but also for CC (due to the 
prohibitive cost of such calculations) as well as for the full- 
classical methods consists of providing pressure broadening 
coefficients without performing the thermal average over 
the distribution of relative kinetic energy of collisions. 
Within this approximation we consider that cross-sections 
are independent of translational energy and approximate 
the partial pressure broadening coefficient y(j^,j 2 ,T) by its 
value at the mean relative speed v, leading to [30]: 

r(it J 2 J) =/ vo- <n) 0'i ,J 2 ; Eki„) (6) 

with E ki „ the kinetic energy associated with v. Omitting the 
thermal average, Eq. (4) now reads 


r(j^T)=fv 

1 


> Sj B evenPj B ° * (/a Jb> ^ kin ) 


+ 3 Zj B oddPj B a(n (/yiJB.^kin) 


(7) 


In the following paragraphs, calculated values obtained 
through the use of Eq. (5) will be denoted as ThAv, while 
values obtained with Eq. (7) will be denoted as MThV for 
mean thermal velocity. 

2.3. Quantum dynamical methods 

Both the close-coupling [31-33] and coupled-states 
(centrifugal sudden approximation) [34] methods were 
used. The Appendix A provides the link between the 
diffusion scattering matrix and the pressure broadening 
cross-sections. We also refer the reader to Ref. [23] for 
technical details about dynamical calculations which are 
very similar to the one used in the present study. 

Because of the relatively large reduced mass of the 
C 2 H 2 -N 2 system, the small rotational constants ( B ^ = 
1.1766 cm -1 , B 2 = 1.998 cm -1 ), in particular that of acet- 
ylene, and the large number of terms of the PES, quantum 
dynamical calculations turned out to be very time con- 
suming 1 and demanded a large amount of memory. Con- 
sequently, we were unable to perform CC calculations for 
total energies (E kin + E rot (jf + E rot (j 2 )) greater than 
350 cm -1 , and CS calculations were performed up to only 
700 cm -1 . It is worth comparing this situation with our 
recent studies of C 2 H 2 -H 2 [21] and N 2 -H 2 [22], where the 
same problems of computational capacity were present 
but -due to the nature of these systems - to a lesser extent 
and higher energies could be reached. Also, simple extra- 
polation functions were used in these works when even 
higher energies needed to be accessed. 

In order to compute the pressure broadening coeffi- 
cients of the R(0) line we have only computed the partial 
pressure cross-sections <7 (1) (/ 1 =0,j 2 ,E kirl ) for j 2 <8 and in 
Eq. (4) used averaged values determined from the first 8 
ones to take into account higher populated j 2 states. Fig. 1 
shows these partial pressure broadening cross-sections for 
j 2 < 8 and illustrates the defects of the coupled states 
approximation (CSA); this problem is well known, see e. 
g. [35-38], Thus, and although our quantum dynamical 
pressure broadening coefficients cannot be qualified as 
benchmarks for other methods, the comparison with 
selected experimental values shown in Table 1 illustrates 
that our calculation indeed provides good values for these 
coefficients. We also note that the use of Eq. (7) instead of 
Eq. (4) does not lead to significant changes. 

The same S-matrix files, used to calculate the partial 
pressure broadening cross-sections for the R(0) lines, were 
used to calculate cross-sections for Q lines a m (jr,j 2 ,E kin ). 
For instance, for the = 5,j 2 =0,E kin ) the thermal 

averaging was performed over 62 CC values between 0 
and 310 cm -1 and 9 CS values between 342 and 667 cm -1 
and an extrapolation was made up to about 2500 cm -1 . As 
in [30] for N 2 -N 2 , one may define an ortho- and a para- 
pressure broadening cross-sections. As already observed 
[29,30] these two contributions are very close at 150 and 


1 For instance, at a kinetic energy of 250 cm -1 forj' 2 =2 and j,=0, the 
duration of a MOLSCAT run is about 10 days with 8 cores working in 
parallel on two quad-core processors (Intel Xeon X5550). Recall that 2 
runs are necessary to obtain the right S-matrix elements which con- 
tribute to the partial pressure broadening cross section of the R(0) line. 


20 


F. Thibault et al. / Journal of Quantitative Spectroscopy & Radiative Transfer 142 (2014) 17-24 



Fig-1- Partial pressure broadening cross-sections = 0,j 2 , F. kln ). in A 2 , 
for the R(0) line and j 2 =0 to 7 as a function of the relative kinetic energy 
in cm -1 ; CC values, small symbols, CS values, big symbols. Rotational 
quantum numbers j i and j 2 correspond to C 2 H 2 and N 2 molecules 
respectively. 

298 K. This may be useful for subsequent authors to avoid 
unnecessary calculations. 

2.4. Semi-classical Robert-Bonamy method 

Dynamical calculations were performed, at the mean 
thermal speed associated to the room temperature, using 
the exact trajectory model determined only by the iso- 
tropic part of the PES and potential terms given by Eq. (3) 
up to Li=L 2 = 4 (see for instance Refs. [11,18-20,22] and 
references therein). We have checked that higher multi- 
polar components of the PES do not provide any significant 
contributions to such calculations of the pressure broad- 
ening coefficients. 

The way to determine pressure broadening cross- 
sections from the standard Robert-Bonamy formalism 
[12], whatever the type of transition operators, is well 
established. The method employed here is very similar to 
the one used in Ref. [18] for N 2 -N 2 and therefore will not 
be repeated again. 

Ma and coauthors [20] have refined the above men- 
tioned formalism by considering effects from the line 
coupling for isotropic Raman Q-lines of a linear molecule. 
More specifically, by correctly evaluating the cumulant 
expansion of the Liouville scattering operator (S} = 
exp( - iS i -S 2 ) up to the second order, they are able to 
take into account the line coupling for isotropic Raman Q- 
lines whose initial and final quantum numbers are iden- 
tical. The method developed in [20] cannot be straightfor- 
wardly adapted for R- or P-infrared lines because for such 
lines the above condition is not valid. As a result, one has 
to develop a new method to evaluate the matrix elements 


of the operator exp]^ +S 2 ). Since this subject merits a 
separate derivation that will not be detailed here. This new 
method and its application to the present system are 
presented somewhere else [26], Finally, let us simply 
mention that in the present study, we have not considered 
the imaginary part of (LSt +S 2 ). Because the interaction PES 
has no vibrational dependence, the Sj term becomes zero. 
Given the fact that iS} term is the main component of Im 
(iSi+S 2 ) and Im(S 2 ) is much smaller than Re(S 2 ), [4,8,26], 
we have decided to neglect lm(S 2 ) in our calculations. 

2.5. Full classical method 

The half-widths of isotropic Raman scattering and IR 
(electric dipole absorption) spectral lines of C 2 H 2 per- 
turbed by N 2 were also calculated in the frame of the 
classical impact approach first proposed by Gordon [16,17], 
Translational and rotational motions of C 2 H 2 and N 2 
molecules were modeled using truly exact three- 
dimensional classical (C3D) equations describing the colli- 
sion of two rigid linear molecules (17 first-order Hamilton 
differential equations in body-fixed coordinates [39]). 
Hence, the trajectories are driven by the isotropic and 
the anisotropic parts of the PES. Spin statistical weights of 
perturbing N 2 molecules were taken into account in the 
usual way. The initial intermolecular distance was set to a 
large enough figure (15 A) to exclude starting interaction 
between molecules. Other collision parameters (initial 
conditions for trajectory calculations) were selected via a 
Monte Carlo procedure. Uniform sampling was applied to 
the directions of rotation and initial orientations of both 
molecules in three-dimensional space. The statistical error 
of calculated linewidths (i.e., RMS error of Monte Carlo 
averaging over initial conditions of collisions) was kept for 
different values of temperature T and rotational quantum 
number j generally at a level less than 1%. The impact 
parameter was bracketed within the range 0-12 A. To 
verify that 12 A was a large enough upper limit, we 
conducted additional test calculations at T=150I< with 
somewhat larger values of the maximum impact para- 
meter, namely b max = 13 A and 15 A. They gave nearly 
identical results to b max =\ 2 A, confirming the fact that 
collisions with fa > 12 A have little influence on line broad- 
ening in our case. 

The majority of calculations were performed only at the 
mean thermal velocity v, which greatly reduces CPU time 
(see Eq. (6)). Several tests were conducted confirming the 
validity of MThV approximation in the cases considered. At 
150 K, for instance, for both R and Q lines, relative 
differences between pressure broadening coefficients 
obtained using the MThV approximation with thermally 
averaged values are no more than about 4%. This percen- 
tage diminishes as the temperature increases (see Fig. 2). 
This behavior is confirmed by quantum dynamical calcula- 
tions and allows the use of Eq. (7) for the SC calculations. 

As a second application of these calculations, Fig. 2 
shows that pressure broadening coefficients for isotropic 
Raman spectral lines become closer and closer to IR 
electric dipole absorption lines as j, and T increases. 
This fact is well known (see for instance Refs. [21,29] and 
references therein). 
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h 

Fig. 2. Comparison of thermally averaged (ThAv) and mean thermal 
velocity (MThV) full classical values at 150 K and 373 K for electric dipole 
absorption (EDA) R-lines and isotropic Raman spectra (IRS) Q-lines and 
comparison of pressure broadening coefficients for such types of lines. 

3. Experimental methods 

3.3. Experimental set-up 

In order to carry out an experimental determination of 
the collisional broadening coefficients of the C 2 H 2 -N 2 
system from isotropic Raman data, we recorded high 
resolution Raman spectra of mixtures of these two gases 
at different pressures and temperatures in the region of 
the Q-branch of the v 2 vibrational band of C 2 H 2 , at 
~1975 cm~\ The experimental setup used for the obten- 
tion of these spectra, built around a stimulated Raman 
spectrometer, has been described in detail in prior articles 
(see for example Ref. [29], which reports a similar colli- 
sional broadening study in C 2 H 2 perturbed by H 2 , and 
references therein). While a number of improvements 
have been introduced in the experimental system in recent 
years [40], its essential characteristics remain the same, so 
in this article we will omit a detailed description and 
direct instead the interested reader to the references 
given above. 

As in Ref. [29], the Ar + laser of the spectrometer was 
operated in its 529 nm line, and the frequency of the cw 
dye laser was scanned between approximately 16,935 and 
16,937 cm -1 . This required the use of Rhodamine 6G in 
the dye laser and Kiton Red in the pulsed dye amplifier. 
The frequency difference of the two laser sources thus 
covered the range between ~ 1972.5 and 1974.5 cm -1 , 
where the Q(l) to Q(15) lines of the v 2 rovibrational band 
of C 2 H 2 lie in the Raman spectrum. Q-branch lines corre- 
sponding to higher quantum numbers than j= 15 could be 
observed, but their diminishing S/N ratio (especially in the 
low-temperature spectra) and the fact that they overlap 
with lines from the v 2 +v 5 -v 5 hot band (clearly visible in 
the room-temperature spectra) made it advisable to 
include only lines up to j=15 in the analysis. The output 
power of the cw Ar + laser was set to 500 mW, and the 
peak energy of the amplified dye pulses to ~15 mj/pulse 
at a repetition rate of 30 Hz, corresponding to the repeti- 
tion rate of the Nd:YAG laser used to pump the dye 


amplifier. This new laser is one of the improvements 
introduced in the system in recent years [40], 

Diluted mixtures of C 2 H 2 (3 vol%) in N 2 were used as 
samples. This concentration was chosen as a compromise, 
since it is high enough to provide good Raman signals in 
the range of pressures at which the experiments were 
conducted while, at the same time, being low enough to 
prevent self-broadening from having a noticeable contri- 
bution to the lineshapes in the recorded C 2 H 2 spectra. The 
mixtures were prepared several days in advance in order 
to insure proper mixing. 

Two sets of measurements were conducted, one of 
them at room temperature (298 K) and the other one at 
150 K. The coolable cell described in [29] was used as 
sample container in both cases. In each one of these sets of 
measurements we followed the same procedure: the cell 
was filled with gas mixture at a certain pressure and the 
Raman spectrum of the Q-branch of the v 2 band of C 2 H 2 
was recorded. Then the pressure was reduced through 
controlled, slow pumping and the spectrum recorded 
again. This procedure was repeated several times in order 
to obtain Raman spectra of the sample at different pres- 
sures. For the room temperature measurements these 
pressures were 70, 60, 50, 40, 30 and 20 mbar, while for 
the ones conducted at 150 K they were 35, 30, 25, 20, 15 
and 10 mbar. All the measurements were carried out twice 
(i.e., two independent spectra were recorded at each 
temperature and pressure) in order to average the final 
results. 


3.2. Data reduction 

The spectra thus obtained were analyzed by fitting the 
spectral lines to Voigt profiles. The range of pressures in 
which the experiments were conducted might suggest 
that more complex profiles that include the Dicke narrow- 
ing effect should have been used, since such narrowing 
should already be visible in this pressure range. Indeed, we 
initially performed fits using Galatry and Rautian profiles. 
However, the S/N ratio of the spectra (~30 for the most 
favorable cases) was not high enough for a clear observa- 
tion, let alone quantification of the narrowing effect, and 
all our attempts failed to provide a reliable value of the 
narrowing parameter different from zero. Thus, a Voigt 
profile was finally deemed adequate for our analysis. Any 
possible underestimation of the collisional line widths 
derived from disregarding the narrowing effect should be 
well within their uncertainty intervals, since these are 
rather large due to the noise that prevents us from 
observing the narrowing in the spectra. 

Another concern regarding the use of Voigt profiles is 
the possible presence of line mixing in the spectra. How- 
ever, any line mixing effect that may affect the lineshape 
can be safely disregarded for two reasons. First, as shown 
by Bouanich et al. [13,14] in Q-branches of X<-n IR 
transitions for pressures below 0.1 atm line mixing is 
negligible, a condition fulfilled by our experiments. 
Secondly, since adjacent lines belong to ortho- and para- 
species, interferences are only possible between second 
(or higher) neighboring lines. 
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4. Intercomparisons of pressure broadening coefficients 

4.2. Isotropic Raman lines 

Fig. 3 provides a comparison of experimental pressure 
broadening coefficients of Q-lines of C 2 H 2 in N 2 baths at 
150 and 298 I< with calculated values. Only CC/CS and FC 
calculations were performed at 150 K while the three 
types of calculations were performed at 298 K. At 150 K 
the FC values are inside the large errors bars while the 
quantum dynamical calculations seem to underestimate 
the experimental values. For the latter kind of calculations 
the MThV values obtained at the kinetic energy of 
133 cm using Eq. (6), are slightly larger than the 
thermally CC/CS values (not shown). This may be partly 
due to the fact that the latter are more affected by the 
default of the CS approximation but also by the average 
integration. This difference between MThV and ThAv 
values diminishes at room temperature (our CC/CS values 
tend to be more and more CS values as the collision kinetic 
energies involved are larger and, generally speaking, the 
thermal average is less important as the kinetic energies 
and therefore the temperature increase, see Fig. 2). At 
room temperature, both the CC/CS and FC values are in 
quite good agreement with the experimental data. Now, 
considering the semi-classical values provided by the 
traditional RB method, we see that they are largely over- 
estimated by about 20%. If line coupling effects are 
included in the calculations (see RB method modified by 
Ma et al. [20,26]) the results are drastically improved. 



Fig. 3. Comparison of experimental (stars) pressure broadening coeffi- 
cients for isotropic Raman QQ't) lines at 150 and 298 K with calculated 
values: full-classical (FC) method (squares), CC/CS method (diamonds), 
standard semi-classical method of [18] (sSC, hollow circles) and refined 
semi-classical method of [20,26] (rSC, dotted circles). 
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Fig. 4. Comparison between experimental (stars) room temperature 
pressure broadening coefficients for IR lines [1] and calculated values: 
FC (squares), CC/CS (diamonds), standard semi-classical [18] (sSC, hollow 
circles) and refined semi-classical [26] (rSC, dotted circles) methods. 

4.2. IR lines 

Experimental data at room temperature are numerous 
(see [9]). We have chosen to compare first our calculated 
values with the smoothed ones obtained from measure- 
ments in the R- and P-branches of the v 5 band [1] and 
which are reported in Fig. 4. We remind the reader that 
we were unable, for technical and CPU time reasons, to 
calculate pressure broadening coefficients for such lines 
except for the R(0) line. Here again, we note that the full 
classical values are in quite good agreement with the 
experimental values, especially for the highest j values. 
The traditional RB semi-classical method strongly over- 
estimates the pressure broadening coefficients whereas 
the new method of Ma et al. [20,26] gives the results in 
between the FC values (or the experimental ones) and the 
traditional RB values. 

5. Conclusions 

We have provided new experimental pressure broad- 
ening coefficients for isotropic Raman lines of acetylene in 
nitrogen at 150 and 298 K. Based on these results we can 
claim that the quantum dynamical and full classical 
calculations performed on the ab initio PES of Ref. [15] 
lead to a good agreement with these experimental data. 
The pressure broadening coefficients calculated with the 
full classical method for IR R-lines in our rigid rotor 
approximation as compared with various experimental 
data for such lines are also quite good and confirmed 
previous calculations performed in [9[. Unfortunately, for 
such lines even for a simple system like C 2 H 2 -N 2 , with a 
quite large reduced mass and small rotational constant, it 
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is still too difficult to routinely handle computationally 
close-coupling calculations to allow a complete intercom- 
parison. Alternative methods are therefore required. 

This paper ends a series of our articles [21-23] devoted 
to an intercomparison of quantum dynamical, full classical 
and traditional semi-classical (Robert and Bonamy) form- 
alisms. The standard Robert-Bonamy method [12] is 
unable to provide accurate pressure broadening coeffi- 
cients without adjusting the PES. It is well known [18- 
23,26] that the standard Robert-Bonamy formalism suffers 
from several weaknesses. Among them, there are two 
particularly severe ones. The first is to treat the transla- 
tional motion and the internal motion separately. More 
specifically, the exchanges of their translational kinetic and 
internal rotational energies are not considered during a 
collision process and trajectories are assumed to be inde- 
pendent of anisotropic potential models. The second one is 
to neglect effects from the line coupling in calculations. In 
order to improve the semi-classical method such that it 
can serve as a well-qualified alternative to the quantum 
dynamical formalisms for considering any kind of inter- 
acting molecular pairs, one must directly address these 
weaknesses. As shown in the previous study for self- 
broadening coefficients of N 2 lines [20] and in the cur- 
rently going study for the C 2 H 2 molecule [26] whose 
rotational constant is even smaller than that of N 2 , taking 
into account the line coupling significantly improves 
results derived from the traditional semi-classical calcula- 
tions. This implies that the refinement effort goes in the 
right direction. 
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Appendix A 

Since CC or CS computations of pressure broadening 
cross-sections are very scarce [21-23,29,30,32,41-46] for 
the case of two linear molecules, this appendix recalls a 
few formulae. Following Ben-Reuven [27] and later on 
Coombe et al. [31], Green [32] and Monchick et al. [33], the 
internal motion of the active molecule is disentangled 
from the bath. The rotational motion of the perturbing 
molecule “2” is first coupled with the relative kinetic 
angular momentum f to form j and then is coupled to 

it to provide the total angular momentum J (J = j i + j ). 
This choice is more natural when one is interested to 
relaxation phenomena tied to the active molecule. In this 
coupling scheme the CC generalized pressure broadening 


cross-section is given by 

^(jJbj'iJJbh-Ekin ) 

ia n j b )( j a n j b 

1 Jb J Ja J [ Jb J Ja 
~ <j a (jif')]'\S l °(E kin + Ej a +E j 2 )j a (j 2 /f) <] > 
<jb(J 20 jlS ,b (E kin + Ej b +E h H/bOzOj > *] 


x [5,. S-.-fSf 

1 I2J2 jj f 


(Al) 


where j a , j b characterize a spectroscopic transition with 
rank n of the radiation, k is the wavevector of the relative 
motion: E kin = 2/<k 2 . 

To be compared with experimental values, these cross- 
sections must be summed over final perturber levels j 2 
and averaged over initial bath conditions: 

ff (n) 0j ; T) = I j2 />j 2 ff <n, 0'i J 2 ; T), (A2) 

where we have used a shortcut notation (jj =j a , j h =j a + 1 
for R lines and j b =j a for isotropic Raman lines) and 

tf (n) 0'i ,h ; T) = Xf 2 <T (n, (h .h J 2 ; T). ( A3 ) 

The MOLSCAT code [28] does not provide generalized 
cross-sections (Eq. (Al)) except the full diagonal term 
(J 2 =j 2 ). Moreover, the S-matrix elements are given in 
the traditional scattering scheme where j, and j 2 are first 
coupled to give j 12 . Therefore we have written a post- 
processor code starting with S* ( (j 1 j 2 )j 1 2 ^ I O'j j 2 )/i 2 O as given 
by the MOLSCAT code to obtain matrix elements 

[32], 

Such cross-sections calculations are not at all imple- 
mented in the MOLSCAT code for diatomic-diatomic 
(linear-linear) molecules within the CSA. To this end we 
have used S-matrix elements, S^ m '((j ,j 2 )j 1 2 1 (JiJ'i )J't 2 ) . as 
given by the MOLSCAT code and coded Eq. (96) of Gold- 
flam and Kouri [34] in another post-processor program. 
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